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We investigate the spin diffusion and transport in a graphene monolayer on Si02 substrate by 
means of the microscopic kinetic spin Bloch equation approach. The substrate causes a strong 
Rashba spin-orbit coupling field ~ 0.15 meV, which might be accounted for by the impurities 
initially present in the substrate or even the substrate-induced structure distortion. By surface 
chemical doping with An atoms, this Rashba spin-orbit coupling is further strengthened as the 
adatoms can distort the graphene lattice from sp^ to sp"^ bonding structure. By fitting the Au doping 
dependence of spin relaxation from Pi et al. [Phys. Rev. Lett. 104, 187201 (2010)], the Rashba 
spin-orbit coupling coefficient is found to increase approximately linearly from 0.15 to 0.23 meV 
with the increase of Au density. With this strong spin-orbit coupling, the spin diffusion or transport 
length is comparable with the experimental values. In the strong scattering limit (dominated by 
the electron- impurity scattering in our study), the spin diffusion is uniquely determined by the 
Rashba spin-orbit coupling strength and insensitive to the temperature, electron density as well as 
scattering. With the presence of an electric field along the spin injection direction, the spin transport 
length can be modulated by either the electric field or the electron density. It is shown that the spin 
diffusion and transport show an anisotropy with respect to the polarization direction of injected 
spins. The spin diffusion or transport lengths with the injected spins polarized in the plane defined 
by the spin-injection direction and the direction perpendicular to the graphene are identical, but 
longer than that with the injected spins polarized vertical to this plane. This anisotropy differs from 
the one given by the two-component drift-diffusion model, which indicates equal spin diffusion or 
transport lengths when the injected spins are polarized in the graphene plane and relatively shorter 
lengths when the injected spins are polarized perpendicular to the graphene plane. 

PACS numbers: 72.25.Rb, 75.40.Gb, 72.80.Vp, 71.70.Ej 



I. INTRODUCTION 

Graphene is considered to be a promising candidate 
for the spintronic apphcations recently, ^"^^ partly due 
to the perfect two dimensionality, gate-voltage-tunable 
charge carrier type and density,'^'* high mobility^"^'^'^'^* 
as well as the potentially long spin relaxation time 
limited by the small intrinsic spin-orbit and hyperfine 
couplings. From the high mobility and long spin 
relaxation time, a long spin relaxation length, favor- 
able to the spin information transport and manipula- 
tion, is anticipated. However, both the spin relaxation 
time and transport length were experimentally found to 
be much smaller than expected. ^'^^'^^"'^^ This suggests 
that the spin relaxation in the experiments is most likely 
to be contributed by extrinsic factors such as the possi- 
ble impurity-enhanced spin relaxation^^'^'' via the EUiot- 
Yafet'^'^ mechanism or the enhanced Rashba spin-orbit 
coupling field^'^^ from the impurities. ^^''^^''^^ The former 
case may exist in a highly dirty graphene sample and 
causes the spin relaxation time Ts to be proportional 
to the momentum relaxation time Tp.^*''^" However, for 
the latter case, the Dyakonov-Perel (DP) spin relaxation 
mechanism'^^ dominates and the relation oc Tp is ab- 
sent. In fact, recently Pi et aJ. reported that Ts increases 
with decreasing Tp in the surface chemical doping exper- 
iment with Au atoms on graphene, '^^ indicating that the 
DP spin relaxation mechanism is important there. How- 



ever, the relation Ts on 1/rp, valid when the DP spin 
relaxation mechanism is dominant and the scattering is 
strong enough, is not obeyed in their experiment.'^^ Nev- 
ertheless, we will show that this deviation can be un- 
derstood by taking account of the strengthening of the 
Rashba spin-orbit coupling with the increasing coverage 
of Au adatoms. The Rashba spin-orbit coupling, referred 
to as an extrinsic one, is due to the broken of the inversion 
symmetry which can be caused by either a perpendicu- 
lar electric field, the interaction with substrate, or the 
atoms adsorbed on the surface. ^'^^'^^''^^''^^ The contribu- 
tion of the electric field to the Rashba spin-orbit coupling 
is small (~ /j,eV under a perpendicular electric field as 
large as 1 V/nm),^^'^^ while the adatoms can effectively 
enhance the Rashba spin-orbit coupling to be of order of 
10 meV by distorting the graphene lattice from sp^ to 
sp^ bonding structure. ^^'^^■'^^ 

In this work, we investigate the spin diffusion and 
transport limited by the DP mechanism in a graphene 
monolayer on Si02 substrate as presented by Pi et al..^^ 
To account for the short spin relaxation time (~70 ps) 
before Au doping in the experiment, '^^ we assume that 
the impurities inevitably present in the substrate, as well 
as the other effects such as the substrate-induced struc- 
ture distortion, cause a strong Rashba spin-orbit cou- 
pling. When the surface chemical doping by Au atoms'^ ^ 
is performed, the Rashba spin-orbit coupling coefficient 
q:r is further strengthened. By fitting the chemical dop- 
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ing dependence of spin relaxation time from the exper- 
imental data,^"'^ we obtain the chemical doping depen- 
dence of aR. It is found that q;r increases approximately 
linearly with the density of adatoms when the latter is 
not too high. With this essential information obtained, 
we then study the spin diffusion and transport in the 
graphene layer. The method utilized in our study is the 
kinetic spin Bloch equation (KSBE) approach which has 
been successfully applied to the study of spin dynamics 
in semiconductors.^^ In the framework of this approach, 
the spatial spin precession frequency during the steady- 
state scattering-free spin diffusion (assumed to be along 
the X-axis) is^''^^^ 

= (2nk + g/XBB)/afe^ek. (1) 

Here fik is the DP term, B is the external magnetic field 
and Ek is the electron energy spectrum. The momentum 
dependence of Wk leads to the inhomogeneous broadening 
in spin precession, with which any scattering (including 
the Coulomb scattering) can cause an irreversible spin 
relaxation along with spin diffusion and transport. 
It is noted that different DP terms as well as different 
energy spectra lead to distinct momentum dependences 
of u;k- For graphene, ek = hvpk with vp = 10^ m/s being 
the Fermi velocity and 

f^k = aR(-sin6'k,cos6'k,0) (2) 

with ^k being the polar angle of momentum k. Therefore 
in the absence of any external magnetic field 

c^k = 2aR(-tan0k,l,O)/(nvF), (3) 

which depends on the angle 0k rather than the magni- 
tude of k. This indicates that the inhomogeneous broad- 
ening is insensitive to temperature and electron density 
as long as is fixed. Therefore the spin diffusion is only 
possible to be modulated effectively by the scattering. 
However, in this work it is revealed that when the scat- 
tering is strong enough (just as in the graphene layer 
under study), the spin diffusion becomes insensitive to 
the scattering. As a result, the spin diffusion is uniquely 
determined by a^- Moreover, the mean spin precession 
frequency (ci;k) = ^jCOjl^O) shows a strong anisotropy 
which can also lead to the anisotropy of spin diffusion 
with respect to the spin polarization direction. This 
anisotropy is found to be quite different from the widely 
believed one predicted from the two-component drift- 
diffusion model.^^"^^ The discrepancy reveals the inade- 
quacy of the two-component drift-diffusion model, espe- 
cially for the cases with spin precession in spatial domain. 

This paper is organized as follows. In Sec. II, we 
present the model and introduce the KSBEs. In Sec. Ill, 
we first investigate the spin relaxation by fitting the ex- 
perimental data from Pi et al.^^ to obtain essential pa- 
rameters and then study the spin diffusion and transport 
in graphene. Both the analytical and numerical investi- 
gations are performed. By comparing the results from 
the analytical and numerical studies, we find that the 



analytical model depicts the zero-electric-field spin dif- 
fusion perfectly and the nonzcro-clcctric-field spin trans- 
port with a small discrepancy which increases with the 
strength of the electric field. At last we summarize in 
Sec. IV. 



II. MODEL AND KSBES 

The n-doped graphene monolayer under investigation 
is on a Si02 substrate perpendicular to the z-axis. The 
depth of the Si02 substrate is assmned to be a = 300 nm 
and the electric field along the z-axis is = Vg/a with 
Vg being the gate voltage. The spins are injected at a; = 
and diffuse or transport along the a;-axis. The external 
electric field, if applied, is along the x-axis, i.e., E = 
Ex. Under the basis laid out in Refs. 25 and 32, the 
Hamiltonian of electrons can be written as^^ 

H = ^ [(ek - M + eEx)6ss' + ■ ct^s'] c^ks^c^ks' 

fi'k.ss' 

+ -ffint. (4) 

Here /i labels the valley located at K or K' point, cr 
denote the Pauli matrices and c^k.5 (c^^ks^) is the an- 
nihilation (creation) operator of electron in /x valley 
with momentum k (relative to the valley center) and 
spin s (s = ±i). Ai is the intrinsic spin-orbit cou- 
pling constant and — e is the electron charge (e > 0). 
The coefficient in the Rashba term J7k [Eq. (2)] reads 
an = (Ez + rj, with the first term contributed by the 
electric field along the z-axis and the second term by the 
substrate (including the impurities initially present in- 
side) as well as the adatoms from surface chemical dop- 
ing. The coefficient C is 5 x 10~^ meV-nm/V (Refs. 26 
and 32). The Hamiltonian H\-a_t consists of the electron- 
impurity, electron-phonon as well as electron-electron 
Coulomb interactions.^^ We adopt the minimal model 
proposed by Adam and Das Sarma"*^ to depict the 
electron-impurity scattering. Within this model, only 
the intravalley electron-impurity scattering is important 
while the intervalley electron-impurity scattering is neg- 
ligible due to the large momentum transfer from one 
valley to the other and the finite distance between the 
impurity layer and the graphene plane. The intraval- 
ley electron-impurity scattering matrix element reads 
|«7k-k'P = Afi|yk-k'pe-2''l'"-'''l,''8 ^Yiere Ni is the effec- 
tive impurity density, d is the effective distance of impuri- 
ties from the graphene layer"^^ and Vk-k' is the Coulomb 
potential under the random phase approximation.'*^ 
The electron-phonon scattering includes the intraval- 
ley electron-acoustic phonon scattering, ^° the intervalley 
and intravalley electron-optical phonon scattering, as 
well as the intravalley electron-optical surface phonon 
scattering. The electron-electron Coulomb scattering 
includes both the intervalley and intravalley scattering, 
with the screening under random phase approximation 
given in Ref. 49. 
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By using the nonequilibrium Green function method, 
the KSBEs are constructed as"^^ 



dt 



dt 



dif 



coll 



dt 



(5) 



Here Pfj_k{x, t) represent the density matrices of electrons 
with relative momentum k in valley p at position x and 



time t. 



dt 



dri 



eE_2fh^i!E2^ are the driving terms 



from the external electric field (the fluctuation of electron 
density is neglected and thus the total electric field is 
taken to be the external one). The diffusion terms due 
to the spatial gradient are 



dpfj.i^{x,t) 



dt 



dsk dpf_,i^ix,t) 



dif 



hdkr 



'Vf cos 0k 



dx 

dpfj.u{x,t) 
dx 



(6) 



dt 



dpf^^{x,t) 



coh 



and £^ 



are the coherent and scat- 



tering terms, respectively. Their expressions can be 

found in Ref. 25. In the steady-state scattering-free spin 

diffusion, the spatial spin precession frequency, given 

by Eq. (3), is immediately obtained according to the 
KSBEs.39-41 



III. SPIN RELAXATION AND SPIN 
DIFFUSION AND TRANSPORT 

In the following, we first study the spin relaxation 
in graphene by fitting the experimental data from Pi 
et al.^^ to obtain information on impurities (including 
the effective density as well as the distance from the 
graphene layer) and the chemical doping dependence of 
the Rashba spin-orbit coupling coefficient. We then use 
the information to study the spin diffusion and transport 
in graphene, first analytically for the case with strong 
electron-impurity scattering only, and then numerically 
with all the scattering explicitly included. 



fitting, the KSBEs are solved in the time domain un- 
der spatial uniform carried out recently by Zhou 
and Wu in the ultraclean graphene monolayer. (An an- 
alytical study of spin relaxation time in graphene is also 
given in Appendix A.) The diffusion coefficient D given 
by Pi et al. is actually for spin instead of charge, al- 
though it is treated as the charge diffusion coefficient 
in the experiment. '^^ In fact, these two coefficients are 
usually close to each other and Jozsa et al. found this 
most likely to be the case in graphene when the elec- 
tron density is high (~ 3 x 10^^ cm~^)'^° due to the weak 
electron-electron Coulomb scattering. Therefore we fit 
the experimental data with the charge diffusion coeffi- 
cient D = ^2^° fivppe, where p-e is the electron mobility 
obtained under a small in-plane electric field. 




2 4 6 

Au deposition (s) 



A. Spin relaxation time 

We fit the chemical doping dependence of spin re- 
laxation time and diffusion coefficient from Pi et al. 
[Fig. 3(c) in Ref. 31] to establish: (i) the density and 
typical distance from the graphene layer of charged im- 
purities initially present in the substrate and those of 
the chemical doping adatoms; and (ii) the dependence 
of aR on chemical doping. The electron density is 
2.9 X 10^2 cm-2 and the temperature T is 18 K.^i The 
electrons are initially polarized in the x-y plane^^ with 
the polarization Pq assumed to be 0.05. To perform the 



FIG. 1: (Color online) (a) Chain curve with triangles: deposi- 
tion time dependence of calculated diffusion coefficient, with 
the Au density growing linearly with the deposition time with 
a fixed rate of 5x10^^ atom/(cm^-s).'^^ Dashed curve with 
closed squares: deposition time dependence of calculated dif- 
fusion coefficient, with the deposition time dependence of Au 
density given by the solid curve (the scale is on the right- 
hand side of the frame). Crosses: experimental data from 
Pi et al..^^ (b) Dashed curve with closed squares: deposition 
time dependence of calculated spin relaxation time, with the 
deposition time dependence of qr shown by the solid curve 
with open squares (the scale is on the right-hand side of the 
frame). Crosses: experimental data from Pi et al..^^ Inset of 
(b): dependence of qr on Au density. 
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We first make use of the value of D « 0.059 m^/s for 
the case without surface chemical doping given in Rcf. 31 
to explore the information of impurities initially present 
in the substrate. This single value of D is not sufficient 
for us to fix both the effective density and distance from 
the graphcne layer of these impurities. However, these 
details are not essential and we just choose two proper 
parameters, e.g., Ng = 2.1 x 10^^ cm~^ and dg = 0.7 nm, 
to ensure D « 0.059 m^/s. The surface chemical doping 
deposits Au atoms on the graphene surface with a growth 
rate of 5x10^^ atom/(cm^-s).^^ By fitting the deposition 
time (adatom density) dependence of D,^^ the distance 
of adatoms from the graphene layer cIau is obtained to be 
about 0.2 nm. Nevertheless, the fitting does not confirm 
with the experimental data well when the deposition time 
exceeds 4 s [compare the fitting data (chain curve with 
triangles) to the experimental data (crosses) in Fig. 1(a)]. 
This indicates that the effective density of adatoms does 
not increase linearly with time any more when the dop- 
ing has been performed for several seconds. Therefore, 
when the doping time is longer than 4 s, we choose the 
proper density of adatoms to reproduce the experimen- 
tal diffusion coefficient. In Fig. 1(a), the deposition time 
dependence of Au density is plotted by the solid curve 
with open squares (the scale is on the right-hand side of 
the frame) and that of the calculated diffusion coefficient 
is shown by the dashed curve with closed squares. 



With the parameters for two kinds of impurities ob- 
tained, we then fit the spin relaxation time to obtain 
q;r under different deposition times. In Fig. 1(b), the 
deposition time dependence of fitted Q!r is shown by the 
solid curve with open squares (the scale is on the right- 
hand side of the frame) and that of the calculated spin 
relaxation time is shown by the dashed curve with closed 
squares. The crosses represent the experimental spin re- 
laxation times under different deposition times. In the 
inset of Fig. 1(b), we also plot the dependence of an on 
Au density Nau- It is shown that qr, increases approxi- 
mately linearly with Au density when the latter is not so 
high. The fitted value of is comparable to the value 
estimated by Ertler et al. when taking account of the 
adatoms, i.e., 0.3 meV."^^ It is noted that aRTp/H has the 
largest value 0.027<1 with Tp y^Neir-J^fie = 0.12 ps 
(Ref. 54) when A'au = 0. Therefore, the electron system 
is in the strong scattering limit (the electron-impurity 
scattering is dominant), let alone when the tempera- 
ture is increased or the chemical doping is performed. 
It is necessary to point out that in the experiment the 
gate voltage Vg is adjusted to keep Nf. constant during 
chemical doping as adatoms also donate electrons to the 
graphene layer. However, Vg does not exceed 200 V and 
the term (E^ = (Vg /a is at least two orders of magnitude 
smaller than aR. Therefore aR » 77 and is solely deter- 
mined by the impurities. When A'^au = 0, 77 = 0.153 meV 
and is contributed by the impurities in the substrate. 



B. Spin diffusion and transport: analytical study 

1. Spin diffusion 

In this section we study the spin diffusion in graphene 
analytically for the case with only the electron-impurity 
scattering. No external electric field is present. We first 
perform the Fourier transformation of the steady-state 
KSBEs with respect to the polar angle and then re- 
tain the equations involving the lowest three orders. ^"'^ 
The neglect of the higher orders will not lose much infor- 
mation in the strong scattering limit where the electron 
distribution approaches isotropy in the momentum space. 
As a result the following second-order differential equa- 
tion about pl^ix) [p^kix) = ^ Jo'' dei,Pi^k{x)e-'^^^ and 
Pnk{x) = Pnk{x,+oo)] is obtained: 



^^kx,[aa.,p°fe(a:;)]] 



(7) 



It is noted that with only the lowest three orders of pj^^ (x) 
considered from the beginning, the electron-impurity 
scattering time is actually absent from the above equa- 
tion (refer to Appendix B for detail). This indicates 
that in the strong scattering limit the spin diffusion be- 
comes insensitive to scattering in graphene. We define 
the "spin vector" as S"j,(a;) = Tr[/9^(x)<T] and S°j^(a;) can 
be solved from Eq. (7) with boundary conditions (refer 
to Appendix B for detail). Then one can calculate the 
total spin signal contributed by all the different electron 
states in two valleys as 



-| /' + 00 p2TT 

S{x) = ^Y.J *y d0j,kTv[pk{x)a 



= - / dkkS%{x). 
Jo 



(8) 



In the following we present the solutions of S{x) under 
three typical boundary conditions. 

For boundary condition (I) S|^fe(0) = (5°fe(0), 0, 0) and 
S°^(-|-oo) = 0, which corresponds to the case with the 
injected spins polarized along the a;-axis. 



S(a;) = 5(0)e- 



Vl + A2 sin( UJX - 


ci sin(a;a;) 



(9) 



For boundary condition (II) S°fe(0) = (0, ^^^(O), 0) and 
S°fc(+oo) = 0, 



S(a;) = 5(0)e- 



X/Iy 



(10) 



5 



For boundary condition (III) S°j.(0) = (0, 0, ^^^.(O)) and 
S°;.(+oo)=0, 



S{x) = 5(0)6 



C2 sin(a;a;) 

^ |. (11) 

-\/l + sin(a;a; — (j)) 



In the above equations 



and 



S{0) = - f dkkSl^iO) 
Jo 

a/7 



{2V2 - 1)^1 + 2^2 ' 



hVF 

2^' 



C\ = - 
8\/2-ll 



(l + V2)Vl+2vl' 



C2 



(20\/2-24)\/l+2v/2 



(12) 

(13) 
(14) 
(15) 
A = 



V7 



and ^ = arctan ^ . It is noted that the spin 



precession frequency given by the simpUfied model is 
w « 1-96^, a Uttie smaller than |(u;k)| = ^ due to 
the approximations made here. 

From Eqs. (9)-(ll) one notices that in the strong scat- 
tering limit, the spin diffusion is not only insensitive to 
the scattering, but also unrelated to temperature T and 
electron density iVg. Nevertheless, the coefficient may 
depend on T and/or A^e, with the relation unclear so 
far. For simplicity we assume q;r, to be independent of 
T and N^. in this work. As a result, the spin diffusion in 
the strong scattering limit is uniquely determined by ccr, 
which is only modulated by chemical doping. Eqs. (9)- 
(11) indicate a strong anisotropy of spin diffusion with 
respect to the spin-polarization direction. For the cases 
with the injected spins polarized along the x- and 2;-axis, 
both the spin signals show an exponential decay in the 
magnitude accompanying with the precession in the x- 
z plane. The spin precessions have the same frequency 
w except for a phase difference. However, when the in- 
jected spins are polarized along the y-axis, the spin sig- 
nal decays exponentially without any precession, i.e., it 
is bound along the y-axis. The above phenomena are 
understood by noticing that the mean effective magnetic 
field felt by the diffusing electrons is along the y-axis as 
(wk) = ^^(O'l'O)- tli^ non-local spin valve experi- 
ments, the spin diffusion length is usually measured from 
the exponential decay of spin signal with the increasing 
spacing between the central spin-injector and -detector 
ferromagnetic electrodes. ^'''^^ In these experiments, the 
ferromagnetic electrodes happen to be magnetized along 
the y-axis and therefore the injected and detected spin 
polarizations are both along the y-axis. With such con- 
figuration, the exponential decay of spin signal with in- 
creasing spacing between the electrodes can be well ob- 
served. However, if the injected spins are polarized in 



the x-z plane, the spatial spin precession is expected to 

be detected. 

Besides the anisotropy of spin precession, the spin dif- 
fusion length also shows an anisotropy as 



It. — ^2 



1.48L 



(16) 



with ly = hv-p/{2aYC)- In fact, when the injected spins 

arc polarized along any other direction in the x-z plane, 
the spin diffusion length is all the same as {Iz) [for this 
case the solution of S(a;) is the combination of Eqs. (9) 
and (11)]. However, based on the widely utilized two- 
component drift-diffusion model'*^"'*^ which gives Ig = 
\/Dts [Eq. (Dll) in Appendix D], one may expect that 
the spin diffusion lengths satisfy 



lx = ly = ^21^ = hVF/{2a-[i) 



(17) 



as the spin relaxation times in time domain follow (refer 
to Appendix A) 



Tx=Ty= 2Tz = h^/{2a%Tp 



(18) 



and D = VpTp/2. It is noted that only when the in- 
jected spins are polarized along the t/-axis, for which no 
spin precession exists, the two-component drift-diffusion 
model gives the result in consistence with that from the 
KSBEs, i.e.. 



hvF/{2an)- 



(19) 



The discrepancy in the anisotropics given by the KSBEs 
and the two-component drift-diffusion model strongly 
indicates the inadequacy of the two-component drift- 
diffusion model. Due to the loss of the off-diagonal spin 
components, i.e., the spin coherence, the two-component 
drift-diffusion model not only fails to predict the spin 
precession in spatial domain in the absence of an ex- 
ternal magnetic field, but also incorrectly inherits the 
anisotropy from the spin relaxation in time domain. We 
emphasize that the reason for the different anisotropic 
properties of spin diffusion in spatial domain and spin 
relaxation in time domain is that the inhomogeneous 
broadening is quite different in these two cases. In spa- 
tial domain the inhomogeneous broadening governing the 
spin diffusion arises from the k dependence of u>k, while 
in time domain from that of ilk- Popinciuc et al. re- 
ported the relationship between the in-plane and out-of- 
plane spin relaxation times directly from the anisotropy 
of spin diffusion via the two-component drift-diffusion 
model. However, based on the above discussion, one 
may realize that studying the anisotropy of spin relax- 
ation time in such a way can be incorrect. Finally, 
from another point of view, if the two-component drift- 
diffusion model is still used, then in order to reflect the 
correct anisotropy of spin diffusion, the spin diffusion co- 
efficient has to differ from the charge diffusion coefficient 
and shows an anisotropy as = 0.5Dz ~ 2.2Dy with 
Dy = v'^Tp/2. 
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It should be pointed out that all the above analy- 
sis and conclusion also apply to the electron system 
where the energy spectrum is parabolic in momentum 
and the linear Rashba spin-orbit coupling term f2k oc 
sin^k, cos^k, 0) is dominant, such as that in the 
asymnicitric InAs quantum wells. That is because the 
steady-state scattering-free spatial spin precession fre- 
quency Wk in this system has the similar momentum de- 
pendence as shown in Eq. (3).'*^ However, for electron 
system in the absence of the DP term but under a mag- 
netic field perpendicular to both the spin polarization 
and spin transport directions such as in bulk silicon^^ and 
symmetric silicon quantum wells, or with the Dressel- 
haus term^^ containing the cubic dependence on momen- 
tum such as in GaAs quantum wells, the situation is 
quite different as u^k depends on the magnitude of mo- 
mentum. In fact, it has been revealed in the symmetric 
silicon quantum wells under an in-plane magnetic field 
that the scattering can suppress spin diffusion effectively 
in the strong scattering limit. 



2. Spin transport 



We further take account of the electric field along the 
X-axis to study the spin transport. Still only the strong 
electron-impurity scattering is included. The second- 
order differential equation about p'^^ix), corresponding 
to Eq. (7) but including the driving term, reads (refer to 
Appendix C) 



hVF 



(20) 



It should be pointed out that when the electric field is 
so large that the electron density matrices p^]^{x,+oo) 
become strongly anisotropic due to the driving of the 
electric field, retaining only the lowest three orders of 
pjife (^) to obtain the above equation of p'^f. (x) may not bo 
sufficient. The second-order differential equation about 
^pfc(^) obtained from the above equation and that 
about S(x) can be obtained by further summing over k 
and p (refer to Appendix C). With the same three differ- 
ent typical boundary conditions presented in the previous 
section, S(x) is solved to have the same form as Eqs. (9)- 
(11) except that the parameters are now electric- field de- 



pendent. Explicitly, 

I' =1' - 

i:.= 



1 



#/2 + F{S) ctR ' 

1 



§12 y4T^74 q;r ' 



/ _ _ 1 Vtr''+48.g^ -1-512 

equations. 



(21) 
(22) 

(23) 
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_ 1 s/g'^+^%g'^+r>Vl 

arctan . In the above 



eE 



In 



1 + 

S'(0)7raR/3fiwF TTe^' 
16V2V<^2 + 7 



(24) 



X Y yjg^ + 48(^2 + 512 - ^2 _^ 8, (25) 
\Jl- <r2/8 + + 48<?2 -I- 512/8. (26) 



Here /5 — l/{kBT) and {p_i) is the chemical potential 
of electrons with spin parallel (antiparallel) to the spin- 
polarization direction. It is noted that when the electric 
field is absent, i.e., (f = 0, all the above solutions recover 
those presented in the previous section. 

In most conditions (such as in the present work) elec- 
trons in graphcne are highly degenerate. In the de- 
generate limit with small spin polarization, S" f« ankip ' 
where fcp = V'^^e is the magnitude of the Fermi mo- 
mentum of unpolarized electrons with density being Ng 
(Appendix C). Differing from the spin diffusion without 
electric field, the spin transport becomes sensitive to elec- 
tron density as S' depends on the electron density. In the 
nondegenerate limit, S' « "^^'^ (Appendix C) and the 
spin transport becomes sensitive to temperature rather 
than electron density. Moreover, with this value of <f , 
Eq. (22) becomes 



e^/3/2 + ^JJ^^/A + lJ}^ , (27) 



where ly is the spin diffusion length without electric 

field [Eq. (14)]. This result recovers that from the two- 
component drift-diffusion model, which apparently fails 
to correctly reflect the anisotropy of spin transport. 
Therefore, our investigation again indicates that only 
when the spatial spin precession is absent, the two- 
component drift-diffusion model gives the appropriate 
depiction of spin transport. 



In Fig. 2 we plot the dependence of l'^ 



and 



on S' . From Fig. 2(a), one notices that the spin trans- 
port length decreases with increasing # (# « ) ■ On 
one hand, this means that when the electron (fensity is 
fixed (e.g., = 10^2 cm^2^ for which the variation of £' 
from —8 to 8 corresponds to a variation of E from about 
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—2.2 to 2.2 kV/cm), the spin transport is suppressed (en- 
hanced) by increasing the electric field parallel (antipar- 
allel) to the spin injection direction. On the other hand, 
this also means that when the non-zero electric field par- 
allel (antiparallel) to the spin injection direction is fixed, 
the spin transport is enhanced (suppressed) by increasing 
electron density. Fig. 2(b) and (c) indicate that the spin 
precession frequency w' and the phase angle cf)' vary with 
S" marginally (with a variation ~ 2 %). In fact, when \S'\ 
becomes even larger, both lo' and 4>' quickly saturate (w' 
approaches and cf)' approaches 7r/2). Therefore, the 
spin precession pattern in spatial domain is insensitive to 
the electric field or the electron density. 



1.98 
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FIG. 2: (Color online) The dependence of (a) spin transport 
length l'x,y,z, (b) spin precession frequency ui' and (c) phase 
angle 0' on S. 



and the two-side injection boundary conditions'^^ are 



-!^n-cr,(31) 



<o — 



^kt 



^ki 



(32) 



Here the injected spins at left boundary a; = are as- 
sumed to be polarized along fi with polarization Pq = 
0.05. X = L stands for the right boundary with L much 



longer than the spin diffusion or transport length. ^ i 
are the Fermi distributions of electrons at the two bound- 
aries when the external electric field is absent. When the 
electric field is present, F^^^^ then stand for the drifted 
Fermi distributions of hot electrons. In the previous 
analytical study the boundary conditions are in fact ap- 
proximated as the single-side injection case. This ap- 
proximation works well when the scattering is strong. '^^ 
By numerically solving the KSBEs, the steady-state dis- 
tribution of spin polarization along fi is obtained as 
Pix) = EpkTrKk(a;,+oo)n • cr]/X;^^Tr[p^k(a::,+oo)] 
and then the spin diffusion or transport length is deter- 
mined from the exponential decay of P{x) (or its enve- 
lope) along the x-axis. 




10" 



T=18 K, N^=2.9xl0'^ cm"^ j 



10 15 
X (|im) 



20 



25 



C. Spin diffusion and transport: numerical results 

The KSBEs need to be solved numerically in order to 
take full account of all the different kinds of scattering as 
well as the large electric field. To numerically solve the 
KSBEs, the initial conditions are set as 



FIG. 3: (Color online) The absolute value of steady-state 
spin polarization \P\ versus position x with the injected spins 
polarized along the x-, y- and z-axis, respectively. The 
squares, circles and triangles are obtained by numerically 
solving the KSBEs with T = 18 K, Ne ^ 2.9 x 10^^ cm-^, 
iVs = 2.1 X 10^^ cm"^ and A'^au = 0. The solid, dashed and 
chain curves are calculated from Eqs. (9)-(ll) with P{x) = 
&(a;)/iVe, Sy{x)/Ne and S^{x)/Ne, respectively. 



P/.k(0,0)- \ n-(T, (28) 

1. Anisotropic spin diffusion 

pL _|_ pL 

PMk(a:>0,0) = -i^^-^, (29) 

^ As revealed by the analytical model, the spin diffusion 

^T^r[pf_,i^{0,0)h ■ cr]/ ^Tr[p^k(0,0)] = Pq, (30) shows anisotropic properties with respect to the polar- 
pk ization direction of injected spins. In Fig. 3, we show the 
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spatial distribution of the absolute value of the steady- 
state spin polarization \P\ for the cases with the injected 
spins polarized along the x-, y- and z-axis, respectively. 
^Au — 0, with which or — 0.153 meV. The squares, 
circles and triangles are obtained by numerically solving 
the KSBEs while the solid, dashed and chain curves are 
calculated by Eqs. (9)-(ll). When or = 0.153 meV the 
analytical model gives Ix = Iz ~ 3.18 ^m, ly w 2.16 /im 
and u! « 0.45 /zm~^. The anisotropy of spin diffusion is 
clearly shown in this figure. It is noted that the simplified 
analytical model almost perfectly recovers the numerical 
results [except that the spin precession frequencies for 
both cases with the injected spins polarized along the x- 
and z-axis are numerically shown to be closer to 
rather than 1.96^^ given by the analytical study (the 
difference is expected from the approximations made in 
the analytical analysis)]. In fact, further numerical cal- 
culations show that varying T from 18 to 300 K and/or 
Ne from 0.5 to 2.9 x 10^^ cm"^ changes the numerical 
results marginally. This is consistent with the conclu- 
sion from the analytical model, i.e., the spin diffusion of 
electrons in graphene is insensitive to T and in the 
strong scattering limit. As a result, in the strong scat- 
tering limit, one can depict the spin diffusion quite well 
with the single parameter ofR via Eqs. (9)-(ll). 

2. Chemical doping dependence of spin diffusion 

In Fig. 4, we plot the deposition time dependence of 
spin diffusion length with n = x, y and z respectively 
by the solid curves. The spin diffusion lengths are di- 
rectly obtained from Eqs. (13)-(14). It is shown that 
with the increase of chemical doping time, or increases 
and the spin diffusion length decreases. For comparison, 
we also plot the deposition time dependence of spin dif- 
fusion length given by the two-component drift-diffusion 
model (chain curves), i.e., Ix — ly — V^lz — \/Dtx, with 
D and Tx given in Fig. 1. The comparison between these 
two sets of results shows that, only when the injected 
spins are polarized along the y-axis, the two-component 
drift-diffusion model yields the same result as that from 
the KSBEs, just as revealed in the analytical study [refer 
to Eq. (19) and the discussion there]. 

3. Effect of scattering on spin diffusion 

The electron system under investigation is always in 
the strong scattering limit and therefore the spin dif- 
fusion becomes insensitive to scattering. However, the 
properties of spin diffusion in the weak scattering limit 
can be different. In order to investigate the spin diffusion 
with scattering strength ranging from the weak to strong 
scattering limit, we artificially vary the impurity density 
in the substrate from to 10"'^^ cm~^. At the same time, 
the chemical doping is absent (no adatom) and aR is kept 
as a constant, e.g., 0.153 meV. We choose T = 50 K, 




Au deposition (s) 

FIG. 4: (Color online) Deposition time dependence of spin 
diffusion length with the injected spins polarized along the 
y- and z-axis, respectively. The results from the KSBEs 
(solid curves) and the two-component drift-diffusion model 
(chain curves) are both plotted for comparison. 

A^e = 5 X 10^^ cm"'^ and n = y. In Fig. 5 we plot the de- 
pendence of spin diffusion length ly on the impurity den- 
sity by the dashed curve. For comparison, we also plot 
the corresponding dependence of spin relaxation time Ty 
on the impurity density by the solid curve (the scale is 
on the right-hand side of the frame). It is seen that with 
the increase in Af^, while ly decreases obviously in the 
weak scattering limit (A^s ^ 0.05 x 10"'^^ cm~^) and then 
saturates in the strong scattering limit, Ty first decreases 
in the weak scattering limit (refer to the inset for detail) 
and then increases almost linearly in the strong scatter- 
ing limit. The two-component drift-diffusion model is 
able to capture the dependence of spin diffusion length 
on Ns by means of the relation ly — yj Dry: while 
D (X Tp (X 1/As, Ty decreases with in the weak scatter- 
ing limit and cx Ns in the strong scattering limit; there- 
fore ly first decreases with Ns and then becomes insen- 
sitive to Ns (the insensitivity of ly to Ns in the strong 
scattering limit is revealed previously by the analytical 
study). It should be emphasized that in Fig. 5 the results 
are shown with aR being a constant. In reality, when one 
further takes account of the increase of q;r with increas- 
ing As, ly should always decrease with increasing A^^, 
from the weak scattering limit to the strong scattering 
limit. 



4- Spin transport under the electric field 

At last we investigate the spin transport under an elec- 
tric field along the x-axis. T = 300 K, A^^ = 10^^ cm'^, 
Ns = 2.1 X 10^2 cm-^ and A^au = 0. The injected 
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FIG. 5: (Color online) The impurity (in the substrate) density 
dependence of spin diffusion length (dashed curve) and spin 
relaxation time (solid curve with the scale on the right-hand 
side of the frame). The inset shows the detail of the solid 
curve in the small density regime. The injected spins are 
polarized along the j/-axis. T = 50 K, A^e = 5 x 10^^ cm^^ 
and qr = 0.153 meV. 
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FIG. 6: (Color online) The absolute value of the steady-state 
spin polarization |P| versus position x under different electric 
fields. The electric field dependence of spin transport length 
Iz is also plotted with the scale on the right-hand side and 
top of the frame, where the squares and double-dotted chain 
curve are obtained from the numerical calculation and from 
Eq. (21), respectively. T = 300 K, = 10^^ cm~^ Ns = 
2.1 X 10^^ cm"^ and A^au = 0. 



IV. CONCLUSION 



spins are polarized along the z-axis. In Fig. 6 the po- 
sition dependence of |P| under different electric fields as 
well as the E dependence of Iz (squares with the scale 
on the right-hand side and top of the frame) are plot- 
ted. It is shown that while the spin-precession pattern 
almost keeps the same with varying E, the spin trans- 
port length is increased (decreased) by increasing the 
electric field along the —x (x)-direction.'^^''''' These re- 
sults are in consistence with the analytical study pre- 
sented in Sec. Ill B 2. For comparison, we further plot the 
E dependence of Iz from Eq. (21) by the double-dotted 
chain curve with the scale also on the right-hand side 
and the top of the frame. It is shown that the analytical 
model depicts the spin transport in the low electric-field 
regime well except when the electric field antiparallel to 
the spin-injection direction is large (e.g, a discrepancy 
reaches 20% when E reaches —2 kV/cm). 

The electron density dependence of spin transport is 
also investigated. In Fig. 7, we plot the density de- 
pendence of spin transport length under the electric 
field parallel {E = 0.3 kV/cm) and antiparallel {E = 
—0.3 kV/cm) to the spin transport direction in (a) and 
(b), respectively. The squares are from the numerical 
calculation and the curves are from Eq. (21). It is clearly 
shown that for the cases with opposite directions of the 
electric field, the density dependences of spin transport 
length have opposite tendencies. 



In conclusion, wc have investigated the spin diffusion 
and transport in graphene monolayer on Si02 substrate 
as presented by Pi et al.,^^ by means of the KSBE ap- 
proach. The substrate (including the impurities initially 
present) contributes a Rashba spin-orbit coupling field 
much stronger than the one modulated by the electric 
field perpendicular to the graphene layer. By surface 
chemical doping with Au adatoms, the Rashba spin-orbit 
coupling coefficient q;r, is increased. By fitting the chem- 
ical doping dependence of diffusion coefficient and spin 
relaxation time,'^^ we obtain the information on impu- 
rities as well as the chemical doping dependence of aR. 
Our fitting finds that q;r increases linearly from 0.15 to 
0.23 mcV with increasing Au density when the latter is 
not so high. With the necessary parameters obtained 
from fitting, we investigate the spin diffusion and trans- 
port in graphene both analytically and numerically. 

The analytical study with only the electron-impurity 
scattering included reveals that in the strong scatter- 
ing limit (just as the situation under investigation in 
the present work), the spin diffusion is uniquely deter- 
mined by ttR. When the injected spins arc polarized 
along the x-, y- and z-axis, the spin diffusion lengths 
are given by the analytical study with an anisotropy as 
Ix = Iz ^ 0.74/iWF/aR and ly = O.Sfivp/aR. Meanwhile, 
the spatial spin precession is present when the injected 
spins are polarized in the x-z plane but absent when the 
injected spins are polarized along the y-axis. Further 
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FIG. 7: (Color online) Electron density dependence of spin 
transport length under electric fields with opposite direc- 
tions: (a) E = 0.3 kV/cm and (b) E = -0.3 kV/cm. The 
squares are from the numerical calculation while the curves 

and 



are from Eq. 

iVAu = 0. 



(21). T = 300 K, Ns = 2.1 x 10^ 



numerical calculations with all the scattering explicitly 
included show that the analytical model depicts the spin 
diffusion pretty well. 

It is noted that the anisotropy of spin diffusion length 
from the KSBEs differs from the one from the two- 
component drift-diffusion model where Ix = ly = V^lz — 
O.bhvp/an. The qualitative discrepancy indicates the in- 
adequacy of the two-component drift-diffusion model due 
to the neglect of the off-diagonal spin components, i.e., 
the spin coherence. In fact, only when the injected spins 
are polarized along the y-axis and the spatial spin preces- 
sion is absent, the two-component drift-diffusion model 
gives the same spin diffusion length as the KSBE ap- 
proach does. 

The analytical and numerical study of spin transport 
under an electric field parallel or antiparallel to the spin 
injection direction is also investigated. In the presence 
of the electric field, the analytical model depicts the spin 
transport with a small discrepancy which increases with 
the strength of the electric field. It is shown that when 
the electric field is applied, the spin precession in spatial 
domain for the cases with the injected spins polarized 
along the x- and 2;-axis remains almost unchanged. How- 
ever, the spin transport length is increased (decreased) 
by increasing the magnitude of the electric field when it 
is antiparallel (parallel) to the spin transport direction. 
Moreover, in the presence of the electric field, the spin 
transport becomes sensitive to the electron density, dif- 
fering from the case of spin diffusion. The spin transport 



is enhanced (suppressed) by increasing electron density 
when the electric field is parallel (antiparallel) to the spin 
injection direction. 
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Appendix A: Spin relaxation in graphene 

We consider spin relaxation in graphene under the spa- 
tial uniform case in the absence of the electric field. We 
only include the electron-impurity scattering. The KS- 
BEs, Eq. (5), are then simplified to be 



i9tp^k(t) 



--[S2k ■ (T,Pf,k(t)\ - — y Mk-k'-?^kk' 

k' 

X (5(ek - EkOb^kW - Ppk'W]- (Al) 



|C^k-k'P 



\U^\,\'^ is the total electron- 



Here Mk_k' 

impurity scattering matrix element contributed by im- 
purities in the substrate and Au adatoms. /kk' = ^[1 + 
cos(0k — fi'k')] is the form factor. By expanding p^-k{t) 
as Pfj.k{t) = y^.i p'..i.(t)e^^^^ . one comes to 



ttR I 



Or, 
2V 



where a± = Cj; ± zu. 



1 k{l - S, 



10) 



and 



(A2) 



d6iMq(l-hcos6')(l-cos;6i) (A3) 
It is noted 



with Mq depending only on |q| = 2fcsin| 
that 4 = 

Retaining the lowest three orders of i-^-, ^ 



0, 



±1, and using the initial conditions pj,;j(0) — (^mP^felO), 
one obtains the second-order differential equation about 
as 



a 



2^2 



M'^yMv.Plkm^^ 



(A4) 



with an affiliated initial condition dtp^^i.{Q) = 0. Defining 
the spin vector as S°j.(i) = Tr[p°^(i)cr], one can obtain 
an equation satisfied by S°^(t) directly from the above 
one, which reads 



dt + -dt + 



(l + '^o.) ^°fcji) = (A5) 
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with a = x,y, z. With the initial condition 9(5°^^ (0) = 
0, S'^^^,^(t) is solved to be 



+ 



e "fc 



-rT-(i+\/i^) 

e 



(A6) 



where = 2-\/2(l + Saz)aRTl/h. When the scattering 
is strong enough and hence <C 1, 



(A7) 



As a result, for spins polarized along the x- and 
y-axis the spin relaxation times are Tx = Ty = 

/{2a\Tl), while for spins polarized along the ^;-axis 
Tz = fi^/(4a|r^). From Eq. (A3) one notices that 
is in fact the momentum relaxation time Tp{k). For the 
highly degenerate electron system in graphene, Tp{k) w 
Tp{k-p) « Tp. Therefore we have Tx = Ty = 2tz = 
h^/{2alTp). 



Appendix B: Spin diffusion in graphene 

The spin diffusion in the absence of an electric field 
is also investigated for the case with only the electron- 
impurity scattering included. Performing angle expan- 
sion on the steady-state KSBEs in a way similar to that 
shown in Appendix A, one arrives at 

io=±l 

2 pUM_„ 



+ 



VY Tl 



(Bl) 



with 7 = an/{hv-p)- Retaining the lowest three orders of 
p''^k(x) one obtains three equations involving pj^^^ix) as 

E P%i^)+^['^+^pUx)] - l[cr-,p;l{x)] = 0,(B2) 

(o=±l 

d.pU^) - lW-,P%{x)] + = 0, (B3) 



9xplk{x)+-f[a+,plk{x)] + 



2 Pji^) 
VF tI 



0. 



(B4) 



From these equations one immediately arrives at Eq. (7) 
with being irrelevant. By multiplying cr and perform- 
ing trace on both sides of Eq. (7), one gets the equation 
satisfied by S°^(a;) which can be written as 

9^-47^ -479, \ (Sl^^{x)\ 

92 _ 4^2 ^o^^(^) ^o_(B5) 

475, 52 -87V KsluA^)) 

With specified boundary conditions, S^j.(x) is solved and 
the total spin signal S(.x) is obtained by Eq. (8), as pre- 
sented in Sec. IllBl. Explicitly, taking the boundary 
condition (I) given in Sec. Ill B 1 as an example, one ob- 
tains S°j.(a;) as 



(Vl + A2 sin(w.T + 4 
I , (B6) 

Ci sin(u;a;) 

with the parameters Ix, w, A, (j) and ci given in 
Sec. Ill B 1. By further summing over k and fi one arrives 
at Eq. (9). 



Appendix C: Spin transport in graphene 

The analytical study of spin transport is carried out 
analogly. The driving term from the electric field in the 
steady state is approximated as 

eE dpf,^{x) _ eE dp^i,{x) de^ 
h dkx h 9ek dkx 

dp%ix) 



eEv-p cos ^k- 



9ek 



(CI) 



Then the Fourier transformation of the steady-state KS- 
BEs reads 

/o=±l 



eE 



-{Si-i + 6n) + 



VF Tl 



0. (C2) 



From this equation one comes to Eq. (20) by retaining 
the lowest three orders of ^^^^(a;). The equation satisfied 
by S%ix) is then 



d^x - eEdxde, - 472 


47^^ - 2eE^d,^ 







dl - eEdxde^ - 472 

dl 



-A^dx + 2eE^de^ 




eEdxd,^ - 87^ 



= 0. 



(C3) 
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Having the experience of solving Eq. (B5), we assume 
that S°^.{x) has the solution as S°^(a;) = S'°^(0)T(a;) 
and therefore S(a;) = S'(0)T(a;). Performing summation 
over /i and k on both sides of the above equation and 
using the trick 



r+oo 

/ dkk[deX,mT{x) 

Jo 



one obtains the equation satisfied by S(a;) as 



S{x), (C4) 



di + 7^a, - -47a, - 272,? 

dl+ 7^a^ - 472 

47a, + 272^ dl+ 7^a^ - 872 



S:c{x) 

X I Sy{x) I = 0, 

S,{x) 



(C5) 



in which S" is given by Eq. (24). With the three differ- 
ent typical boundary conditions presented in Sec. Ill B 1, 
S(a;) is solved to have the same form as Eqs. (9)-(ll) 
except that the parameters are now given in Sec. IIIB 2. 

We now calculate S in both the degenerate and 
nondegenerate limits. In the degenerate limit, 

/3/iWF\/^(\/l + -Po - Vl - Pa) ~ PhypVTrKPo- Mak- 
ing use of the relation S{0) = NgPo, one has 



eE 



eE 



\/Tr NeaB. ctnkF ' 



(C6) 



where fcp = -v/TrA^e, the magnitude of Fermi momentum 
of unpolarized electrons with density being N^. In the 
nondegenerate limit, In ^]||''^^^ e'^^'f — e^'^^ and 



1 /■+°° 
S{0) w - / dkk[i 

Jo 



therefore 



eE^hv-p 



(C7) 



(C8) 



Appendix D: Derivation of two-component 
drift-diffusion equation from KSBEs 

The two-component drift-diffusion equation can be de- 
rived from the KSBEs in the collinear spin space^'^ with 
the z-axis along the initial spin-polarization direction n, 
by neglecting the spin coherence (i.e., the off-diagonal 



components of the density matrices). The density ma- 
trices then have the diagonal form as 5[//ikt(2;, i) + 
fu.H{x,t) + (^kt(a;,t) - fu,ki.{x,t))a^]. In the following 
we present a brief derivation of the two-component drift- 
diffusion equation from the KSBEs with only the strong 
electron-impurity scattering considered. Other kinds of 
scattering can also be incorporated similarly under elas- 
tic scattering approximation. The spin relaxation time 
is Tfi and the momentum relaxation time is Tp (both are 
given in Appendix A) and we neglect their momentum 
dependence hereafter. The two-component drift-diffusion 
equation is obtained from the equation of continuity and 
the equation of current, both to be derived from the KS- 
BEs. 

The equation of continuity is derived as follows. By 
multiplying the KSBEs [Eq. (5)] with ^(liba^) and then 
performing the trace, one obtains the simplified KSBEs 
for each spin band (tr =t, i) as 

dfij,ka{x,t) _ eE df^k„{x,t) ^ ^.Q„g^^ ^//^ka(a;,t) 
dt h dkx dx 



f IJ,k.<7 {x ^ t) fj^]i—a{x^t) 

2^ ■ 



(Dl) 



The right-hand side of the above equation comes from 

the term Tr{i(l ± a-2)[5f/9^k(a;, t)|coh + 5tp^k(a;, t)|scat]}, 
which can be calculated with the aid of the KSBEs in the 
time domain [Eq. (Al)] and the corresponding solution 
[Eq. (A7)]. Performing summation over // and k on both 
sides of Eq. (Dl) in the steady state, one comes to 

2. + ^ 2^^^ cosek^k.(a;) 



N^{x)-N_^{x) 
2ta 



(D2) 



where N^{x) is the electron density with spin a at posi- 
tion X. Up to the first order of the electric field E, the 
first summation over k in the above equation leads to 
zero when /^uko- is approximated by f^i^„, the distribu- 
tion in equilibrium. Defining the charge current along 
the X-axis with spin a as 



Ja{x) = -^evxfij,i,^{x) 



(D3) 



with Vx = vf cos ^k, one has the equation of continuity 
1 dMx) _ N^{x) - N_„{x) 



dx 



2rfl 



(D4) 



We then calculate the current from the diagonal 
part of the KSBEs 

dff,ka{x,t) eE df^\^a{x,t) dff,ka{x,t) 

+ Vf cos t7k " 



dt h dkx 

^ fnkaix,t) - f°i^^{x,t) 
Tp 



dx 



(D5) 



13 



where the right-hand side of the equation conies from the 
electron-impurity scattering (Appendix A) . In the steady 
state, multiplying —evx on both sides of the equation and 
then summing over and k, one comes to 

^ E ^^^^^ - E ^>^^%^ 



fik. 
Ja{x) 



dx 



(D6) 



Again, up to the first order of the electric field, one has 
J^{x) = eiiaEN^{x) + eDd^N^ix), (D7) 



where the mobility /ig 



and the charge diffusion 



coefficient D = ^VpTp. For the case with small spin po- 
larization, m = He = ^yj^ in which = N^-\-N\_ 
is the total electron density. 

Finally, the two-component drift-diffusion equation is 
obtained by combining Eqs. (D4) and (D7) 



N,{x) - N_,{x) 



dx 



dx"^ 



2ta 



(D8) 



This equation is consistent with that in the 
hterature.'*^"'*'^ The equation of AA^ = - 
then reads 



dAN{x) ^d^AN{x) _ AN{x) 



dx 



dx'^ 



(D9) 



With boundary condition AN{+(X)) — 0, AN{x) is 
solved as AN{x) = AA''(0)e~^/'* where the spin trans- 
port length^^ 



l^eE bfJ.eE.2 1 



When the electric field is zero. 



la = y/Dra. 



(DIO) 



(Dll) 



As indicated by Eqs. (DIO) and (Dll), in the frame of 
the two-component drift-diffusion model, the anisotropy 
of the spin transport or diffusion is solely determined by 
that of the spin relaxation and no spatial spin precession 
can be obtained. 
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